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THERAAAL CONDUCTIVITIES OF ORDINARY AND ISOTOPICALLY SUBSTITUTED 


POLAR GASES AND THEIR EQUIMOLAR MIXTURES* 
by Charles E. Baker and Richard S. Brokaw 
Lewis Research Center 

SUMAAARY 

Thermal-conductivity measurements have been carried out over the temperature 
range 300° to 475° K using an eight-cell hot-wire thermal -conductivity apparatus. The 
thermal conductivities of the deuterated gases are greater than those of their unsubsti- 
tuted counterparts at the higher temperatures (presumably because of their higher heat 
capacities). Conductivities of the equimolar mixtures lie about midway between the con- 
ductivities of the pure components. Results are analyzed to test a recent hypothesis that 
a resonant exchange of rotational quanta causes the thermal conductivity of a polar gas to 
seem anomalously low in relation to its viscosity. These results are consistent with this 
hypothesis, the only qualification being that near-resonant exchange is also important. 


INTRODUCTION 

The thermal conductivities of highly polar gases such as hydrogen fluoride, water, 
and ammonia appear to be anomalously low in relation to their viscosities. Mason and 
Monchick (ref. 1) have suggested that this effect is largely a resvilt of a resonant ex- 
change of rotational energy, presumed probable on grazing self-collisions of polar mole- 
cules. Hence, a grazing collision with exchange is equivalent to a head-on collision 
without exchange insofar as the transport of the rotational quantum is concerned. Baker 
and Brokaw (ref. 2) have recently measured the thermal conductivities of gaseous water 
(HgO) and deuterium oxide (DgO), as well as their equimolar mixture, in an effort to 
test experimentally the theory of Mason and Monchick. As has been previously pointed 
out (ref. 2), a polar gas and its deuterated counterpart should be a good test of the theory 
because the substitution of a deuterium atom for a hydrogen atom increases the moment 
of inertia by approximately a factor of 2 but changes the mass only slightly. Since the 
resonant correction depends inversely on the 3/2 power of the moment of inertia, the 
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resonant correction for the deuterated gas will be approximately one-third that of the 
ordinary gas. 

The following conclusions were reached on the basis of the work with the H 2 O-D 2 O 
system (ref. 2): 

(1) The data on the pure compounds seem to confirm the postulate of Mason and 
Monchick regarding the importance of resonant transfer of rotational quanta. 

(2) On the other hand, the data on the mixtures suggest that a resonant exchange of 
rotational quanta is not an important factor in determining the thermal conductivity of 
water vapor. 

(3) It was suggested, tentatively, that exchange of rotational energy on grazing col- 
lisions is important but is of a classical nature. 

These conclusions were not definitive because of inadequacies both in the theory for 
asymmetric top molecules (such as water) and in the theory for polyatmoic gas mixtures. 
In addition, it has been pointed out (refs. 3 and 4) that the H 2 O-D 2 O system was, in fact, 
an equilibrium ternary mixture of H 2 O, HDO, and D 2 O. One way to overcome the first 
shortcoming is to study gas pairs for which the theory does exist, namely, linear dipoles 
such as hydrogen chloride (HCl) and deuterium chloride (DCl), and symmetric tops such 
as ammonia (NHg) and ammonia-dg (NDg). As for the second point, the work by Mason 
and Monchick on the thermal conductivity of polyatomic and polar gases (ref. 1) has re- 
cently been extended to include gas mixtures (ref. 5). 

The thermal-conductivity equations obtained by Mason and Monchick are based on 
the formal kinetic theory of Wang- Chang and Uhlenbeck (ref. 6) and Taxman (ref. 7). 
However, Waldman (ref. 8) has recently pointed out that certain simplifying assumptions 
made in both these theories are not valid; therefore, neither of these theories is actually 
as rigorous as it had previously been considered to be. Consequently, the equations of 
Mason and Monchick must also be regarded as less rigorous. Nonetheless, their formu- 
las for calculating the thermal conductivity of polyatomic and polar gases appear to be 
the best presently available and will be used in the data analysis of this report. 

The experiments reported herein were undertaken to determine whether the tentative 
conclusions reached with H2O-D2O are correct. Conductivities of pure NHg, NDg, HCl, 
DCl, CH^ (methane), CD^ (methane-d^) and equimolar mixtures of HCl-DCl and CH^-CD^ 
were measured in the temperature range 300° to 475° K. The conductivities of 

1C 15 o 

N Hg and of the equimolar mixture of N Hg and NHg were measured at 300. 0 and 

424. 2° K, while the conductivity of the equimolar NHg-NDg mixture was measured at 
300. 0° K only. 
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EXPERIMENTAL 


Apparatus and Procedure 


The hot-wire thermal -conductivity apparatus used for these conductivity measure- 
ments consists of four pairs of carefully matched platinum-iridium filaments mounted in 
a stainless -steel block. The cells are connected as elements of a constant-current 
Wheatstone bridge. When an instrument such as this is used as a thermal-conductivity 
detector in analytical work, as is often the case, a bridge current of 400 to 500 milli- 
amperes is normally used. Under these conditions, the difference between the tempera- 
ture of the hot wires and the temperature of the block, AT, is 300° to 400° C. For this 
work, however, a bridge current of only 100 milliamperes was used, which produces a 
AT ranging from 2° to 25° C, depending on the thermal conductivity of the gas in the 
cells. 

The voltage unbalance produced when gases of unequal thermal conductivity were in- 
troduced into the test and reference cells was measured on a sensitive potentiometer. 
The instrument responds to the reciprocal of the thermal conductivity 


E 




( 1 ) 


where E is the voltage imbalance with a gas of unknown conductivity in the four test 
cells and a reference gas in the other four cells; is the voltage with the reference 

gas in all eight cells; A and Aj,gj. are the thermal conductivities of the unknown and 
reference gases, respectively; and b is a constant characteristic of the apparatus (it 
may be slightly temperature dependent). For this work, nitrogen was used as the refer- 
ence gas, and helium and argon were used as the calibrating gases. 

All the thermal-conductivity measurements were made at atmospheric pressure 
(750 torr). The desired temperatures were maintained by a large oil bath that could be 
controlled to within ±0. 02° C. Temperature measurements were made with a platinum 
resistance thermometer. The gas handling system was conventional in all respects. 
Since all the gases used in these experiments are readily condensed by liquid nitrogen, 
transfer of a gas from a storage vessel to the test cells and back to the storage vessel 
was readily achieved by means of strategically placed cold fingers. 

As previously mentioned, nitrogen was generally used as the reference gas; how- 
ever, all the isotopically substituted gases as well as the equimolar mixtures were mea- 
sured against their nonsubstituted counterparts, that is, NDg against NHg, DCl against 
HCl, etc. Thus, these small electromotive forces were measured with enhanced preci- 
sion. 


3 



Voltages relative to argon are re- 
corded in table I. The measurements 
for each system were not made at 
identical temperatures, since mea- 
surements for a given system were 
made over the entire temperature 
range before changing to another gas 
pair. The experimental voltages were 
reproducible to better than ±0. 005 
millivolt. Because helium and argon 
are monatomic gases, it is possible to 
compute conductivity from viscosity 

through a rigorous, experimentally verified, theoretical relation, to less than a few 
parts per thousand. Thus, thermal conductivities for the calibrating gases were as- 
signed after considering both experimental viscosity and thermal-conductivity data. The 
details of this method have been fully described previously (ref. 2). The conductivities 
selected for calibration are given in table II. 



Figure 1. - Comparison of ammonia data with recent literature 
values. 


Experimental Results 

The thermal conductivities derived from tables 1 and II by means of equation (1) are 
presented in table 111. The table is again subdivided into three parts since the measure- 
ments for the different systems were not made at identical temperatures. Although the 
primary purpose of this work is to measure accurately the relative differences between 
the ordinary gases used and their isotopically substituted counterparts, the absolute 
values obtained are generally in agreement with the literature. This is illustrated in 
figure 1, where the present data on ammonia are compared with the recent measvtre- 
ments of Richter and Sage (ref. 9). The agreement seems satisfactory. 


ANALYSIS AND DISCUSSION 


The expression that Mason and Monchick (ref. 1) derived for the thermal conductivity 
of a polyatomic gas may be written as follows: 


X 


IL\ 

M 



''^trans 



( 2 ) 
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where X is thermal conductivity, rj is viscosity, M is molecular weight, p is density, 

is the average coefficient for the diffusion of internal energy, and 

tx* 3JXS 

are the translational and the internal contributions to the molar heat capacity, is the 
numoer of collisions for relaxation of the internal mode, and Cj^ is the heat capac- 
ity associated with that mode. In small rigid polyatomic molecules, the only collision 
numbers small enough to affect thermal conductivity significantly are those associated 
with rotational relaxation. For nonpolar gases, is commonly assumed to be equal 

to the self-diffusion coefficient, For polar gases, however. Mason and Monchick 

(ref. 1) suggest that 


D 


int ~ 


'11 


1 + 6 


(3) 


where 6 is a correction term calculated from the theory of resonant collisions. Mason 
and Monchick (ref. 1) have developed expressions from which the resonant correction for 
linear dipoles and symmetric tops can be computed. The expression used in the analysis 
of the linear dipoles HCl and DCl is as follows: 


6 = 


<^ 4 >| 


'3 ph V 
16 a^kT 


m 


1/2 



(4) 


Here <a^>, taken to be 0. 44, is the mean value of a dimensionless quantity involving the 
rotational quantum numbers, ji is the dipole moment, h is Planck's constant, is 


( 1 , 1 ) 




1/2 


in the nota- 


the kinetic -theory diameter characteristic of diffusion (Oj^ s a [si 
tion of ref. 10), k is the Boltzmann constant, m is the molecular mass, R is the univer- 
sal gas constant, and I is the molecular moment of inertia. For near spherical top di- 
poles such as NHg and NDg, the expression for 6 is somewhat more complicated than 
equation (4), as a result of the presence of three moments of inertia, two of which are 
equal. This relation is 

. _ Stt / 3 ph Y R % \ 

" 16 OjjkTj i^(i^)l/2 \4 ' 4lJ 


(5) 


where is the moment of inertia about the dipole axis and Ig = are the other two 
moments. For spherical top molecules such as CH^ and CD^, the three moments of 
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inertia are equal. It is easily seen that when = Ig = equation (5) differs from 
equation (4) only by the factor 5n/16, 

The method of analysis of the pure gas data is quite different from the procedure 
used with the equimolar mixtures. Thus, it seems desirable to discuss ^ the pure gas 
results first and then follow this with ^ the mixture results, rather than to discuss 
both pure gas and mixture results for each gas pair before proceeding to the next system. 


Pure Gases 

The experimental data on all the pure gases were analyzed as follows; The neces- 
sary resonant corrections were calculated from equations (4) and (5) and the results were 
used to compute from equation (3), The collision number ^rot was then calculated 

from equation (2) after X had been replaced by the experimentally determined thermal 
conductivity. The ratio of the collision number of the deuterated gas to that of its un- 
deuterated counterpart was then compared with a theoretical value based on the results 
of Sather and Dahler (refs. 11 and 12) for a roi:^h sphere: 


41 


ma 


Z OC 

rot 


IL 
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Here m and I are the same as in equation (4), while a is the viscosity collision 

r (2 ^ 

diameter (a = a ' in the notation of ref. 10) and <j is the diameter of the 

rot^h sphere core. The rough sphere model was chosen because the experimental re- 
sults of O’Neal and Brokaw (ref, 12) for a variety of near -spherical molecules (including 
methane and even ethane) suggest that this model is a reasonable one for substances that 
do not differ markedly from spherical symmetry. 

Some temperature-independent properties necessary for the data analysis are 
listed in table IV. Certain temperature-dependent quantities that are also needed for 
data analysis are shown separately for each gas pair in tables V to VH. The viscosity 

values for ammonia were taken from the best curve through the most recent viscosity 

15 

data available (refs. 13 to 16). Viscosities of NDg and N Hg were calculated from the 
NHg viscosities. Since there were no data available that compared the viscosities of the 
isotopic pairs, the viscosity collision cross sections for all three gases were assumed 
to be the same. Similarly, the viscosities for HCl listed in table VI were taken from the 
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best curve through the data available in the literature (refs. 17 to 20), and the viscosity 
for DCl computed from these data, again with the assumption that the viscosity collision 
cross sections are equal. The viscosities for CH^ were obtained by interpolation from 
the data of De Rocco and Halford (ref. 21). Viscosities of CD^ were computed assuming 
the viscosity collision cross section of CH^ is 2. 2 percent larger than that of CD^. This 
was deduced from the viscosity measurements of van Itterbeek (ref. 22), and suggests 
that the polarizability of CH^ is somewhat larger than that of CD^. A similar effect has 
been observed and discussed for the Hg-Dg system (ref. 23). 

The heat capacity data given in tables V to vn were calculated from spectroscopic 

constants (refs. 24 to 33)^ by a computer program described by McBride, et al. (ref. 

34). The heat capacity values for N^^Hg at 300. 0° and 424. 2° K were obtained from the 

calculated NHg and NDg heat capacities by interpolation on the basis of the square root 

of the reduced mass. As one would expect, the heat capacities for are only 

15 

slightly different from those for NHg; the N Hg heat capacity is larger by just over one 
part per thousand. 

The dimensionless quantity pD^^/rj that involves the self-diffusion coefficient was 
computed as (6/5) (A*) (ref. 10, p. 540). The quantity (A*) is a ratio of collision in- 
tegrals (essentially the ratio of the viscosity and diffusion cross sections) and is quite 
insensitive to temperature and the details of the inter molecular force law, at least for 
spherically symmetric potentials. Since there is no experimental information on the 
self-diffusion of ammonia and only room temperature data for hydrogen chloride, (A*) 
was taken from calculations for a modified (angle independent) Stockmayer potential 
(ref. 35). The values of (A*) used for methane were taken from calculations based on 
an exponential-6 potential (ref. 10, p. 1176)(a! = 14, e/k = 152. 8). The determination of 
(A*) is the only place in the analysis where it is necessary to make any assumption about 
the intermolecular force law. 

All the expressions, molecular properties, temperature-dependent quantities, and 
procedures necessary for the data analysis have now been presented. The first gas pair 
that will be discussed is CH^-CD^. The analysis is the simplest for this system since 
these molecules are spherical tops with no dipole moments; therefore, the resonant cor- 
rection is zero (see eq. (5)). Thus is taken to be pD^^/r] and is computed as 

(6/5) (A*) as before. By substituting the experimentally determined thermal conductivity 
into equation (2), the collision number is easily calculated for both CH^ and CD^. 


^For CH^, refs. 24 to 26; for CD^, refs. 25 and 27; for HCl and DCl, ref. 29; for 
NHg and NDg, refs. 30 to 33. 
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Collision numbers for rotational relaxation calculated by this procedure, as well as 
the ratios of to are presented in table V. The ratio obtained by substitu- 


ing the appropriate quantities from table IV into equation (6) is 



0. 6502. 


The room-temperature collision number for methane is considerably lower than the 
value of 9 collisions obtained by the temperature recovery factor method (ref. 11) which 
also utilizes equation (2). The value of decreases with increasing temperature as 

did the Z^.^^ for H 2 O (ref. 2). It is interesting to note that collision numbers for sulfur 
dioxide (SO 2 ), recently calculated from experimental thermal- conductivity data (ref. 36) 
by means of equation (2), increase with increasing temperature. However, too much 
significance should not be placed on these temperature trends since the collision number 
is extremely sensitive to the absolute value of the thermal conductivity. 

The agreement between the ratios of the experimentally determined collision num- 
bers and the theoretical ratio obtained from equation (6) is quite good at the higher tem- 
peratures, whereas the experimental ratios are somewhat higher than the theoretical 
value at the lower temperatures. It should be emphasized that the ratio of Z numbers 
is extremely sensitive to the ratio of the viscosity cross sections. For example, if one 
assumes that these cross sections are equal (rather than a^jj j<^Qjy = 1.022 as de- 
duced from viscosity data, ref. 21) the ratio of Z numbers becomes about unity. 

In summary, for nonpolar spherical tops, rotational collision numbers obtained from 
experimental thermal conductivity data by means of equation (2) are in reasonable agree- 
ment with current theory for rough spheres, a conclusion previously reached by O'Neal 
and Brokaw (ref. 12). 

As far as polar molecules are concerned, the analysis is simplest for linear dipoles 
such as HCl and DCl. Rotational relaxation collision numbers obtained from the experi- 
mental conductivities and resonant corrections calculated from equation (4) are given in 
table VI. The room temperature value for Z^^^ is in good agreement with the value 
(approximately 7 collisions) recently obtained by Breazeale and Kneser (ref. 37) at 0° C, 
using the acoustic absorption technique. On the other hand, if collision numbers are 
calculated with the assumption that the resonant correction is zero, values are obtained 
(1. 2 to 1. 3 collisions) that are abnormally low when compared with this e3q)erimental 
value. That a collision number of even 6 or 7 for HCl is surprisingly low can readily be 
seen by comparing it with the Z^^^. for nitrogen (N 2 ). The moment of inertia of N 2 is a 
little more than five times greater than the moment of inertia for HCl, whereas their 
masses are about the same. Thus, on the basis of equation (6), one would expect 

Z to be about four times greater than Z . , which is about 5. 5 coUisions 

^°%C1 


(ref. 38). Breazeale and Kneser have attributed the low Z^^^ for HCl to its large dipole 
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moment, and it is probable that the enhanced efficiency with which HCl molecules relax 
rotational energy is due to the long-range dipole-dipole interaction. The collision num- 
bers of DCl are nearly independent of temperature; however, those for HCl clearly de- 
crease with increasing temperature. The agreement between the ratios of the collision 
numbers obtained from experiment and the theoretical ratio is satisfactory. (In con- 
trast, the ratios of the collision numbers calculated with the assumption that the reso- 
nant correction is zero are all unity, clearly in disagreement with all current theories 
on rotational relaxation. 

15 

The final system of pure polar gases to be analyzed includes NHg, NDg, and N Hg. 
Rotational relaxation collision numbers calculated from the experimental data and reso- 
nant corrections calculated from equation (5)^ are given in table vn. The = 2. 3 

collisions obtained for NHg at 300° K is much smaller than the value of 9 collisions 
quoted by Herzfeld and Litovitz (ref. 39) but is in agreement (as it should be) with the 
value of 1. 8 collisions obtained by Mason and Monchick (ref. 1) by fitting thermal con- 
ductivity data. In addition, the NHg value of 2. 3 collisions is consistent with the 4.0 col- 
lisions obtained for CH^ (see table V); their masses and moments of inertia are similar 
but NHg has a large dipole moment (which makes for a small Z^^^), while methane is 
nonpolar. Again, Z^.^^ shows a decrease with increasing temperature, the dependence 
being much the same as that obtained for HgO (ref. 2). The agreement between the ratios 
of the experimentally derived collision numbers and the theoretical ratio is comparable 
to that obtained for the other two systems. Unlike the previous results, however, the 
agreement is better at the lower temperatures than at the higher ones. 

In summary, the data on the pure gases seem to support the conclusion reached on 
the basis of the pure HgO and DgO results; that is, all in all, the data appear to agree 
with the formula of Mason and Monchick for the thermal conductivity of polyatomic gases 
(eq. (2)) and with the correction to this expression necessary for highly polar gases 
(eq. (3)). 


2 

The derivation of equation (5) involves, among other approximations, the replace- 
ment of summations by integrations. In the case of NHg, Mason and Monchick (ref. 1) 
carried out the direct summation on a high-speed computer and found that the analjdic 
expression for the resonant correction gave values that were 5 percent too large at 
273° K and 0. 9 percent too large at 573° K. Thus the ratio of 5 ^eq. (5^ to 6 [direct 
summation] as a function of temperatvtre can be determined from a plot of 
6[eq. (5)]/6 [summation] against h^/(87r^kI^)T. This procedure was applied in comput- 
ing the 6 values (shown in table VII) used in determining Z^^^ from conductivity data. 
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Equimolar Mixtures 


In a recent paper (ref. 2) on H 2 O-D 2 O, the present authors pointed out that Mason 
and Monchick’s postulate regarding resonant exchange of rotational energy could be 
qualitatively tested by comparing the thermal conductivity of two polar gas pairs with that 
of their equimolar mixture. The resonant effect is based on a matching of rotational 
energy levels. Consider, for example, an equimolar mixture of HCl and DCl. In such a 
mixture, one-half the collisions are self-collisions of HCl and DCl in which resonant ex- 
change of rotational quanta may occur. The remaining collisions are between HCl and 
DCl. Since the differences between rotational energy levels in DCl are about half as 
large as the differences in HCl, exchange is not possible, at least for certain collisions. 
Consequently, a mixture might be expected to have a thermal conductivity somewhat 
greater than either of the pure constituents. 

The thermal conductivities of the equimolar CH^-CD^ mixtures illustrate the results 
that would be expected if no resonant effect were involved because CH^ and CD^ are non- 
polar. These results are presented in table VIH along with thermal conductivities of the 
mixtures calculated assuming a linear mixing rule. (The equations of Monchick, Pereira, 
and Mason, ref. 5, were used for a sample calculation at 300° K and agreed with this 
linear average to within one part in ten thousand, ) The experimental conductivities agree 
with the calculated values to within 2 to 3 parts per thousand. The experimental values 
are higher than the calculated values at all temperatures; this suggests that the very 
slight positive deviation from a linear mixing rule may be caused by a slightly smaller 
cross section for the unlike interaction. 

The next system for analysis and discussion is the equimolar HCl-DCl mixture. The 
experimental conductivities will be compared with values calculated using the first-order 
linearized form of the equations of Monchick, Pereira, and Mason (ref. 5) mentioned in 
the INTRODUCTION. This form of their results can be regarded as equivalent to the 
Hirschfelder-Eucken formula for the thermal conductivity of polyatomic gas mixtures 
(refs. 40 and 41) plus a large number of first-order correction terms. 

In calculating the mixture conductivity it is necessary to know the value of tiie reso- 
nant correction for the various interactions in order to compute the coefficient of dif- 
fusion for internal energy (see eq. (3)). For the like interactions, these corrections will 
be the same as for the pure gases; for the vuilike interactions the situation is more com- 
plicated. If the resonance effect actually requires an exact matching of rotational energy 
differences, the correction for the unlike interactions would be zero. There is recent 
evidence, however, that rotational energy exchange does take place in HCl-DCl colli- 
sions (ref. 42), That such an exchange is reasonable can be shown by the following argu- 
ments. The rotational energy of a diatomic molecule is given by 
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(7) 


£rot = W + 1) 


rot 


2 2 

where = h /Stt and J is the rotational quantum number. Since rotational 

transitions are generally restricted to those for which AJ = ±1, the energy difference 
between rotational levels J and J + 1 is 


^J+1 ■ = 



( 8 ) 


The ratio of the moment of inertia of DCl to that of HCl is 1. 9448 (see table IV); there- 
fore, 6 , = 9 , 1. 9448. Thus, the rotational energy level differences for HCl 

and DCl in units of 9 

rotjjci 


jl are as follows: 


HCl 

2 

(J= 1) 


4 

(J = 2) 


(J 

6 

= 3) 

DCl 1. 028 

(J = 1) 

2.057 
(J = 2) 

3.085 
(J = 3) 

4. 114 
(J = 4) 

5. 142 
(J = 5) 

6. 

(J 

170 
= 6) 


etc. 


etc. 


It is evident that, as a consequence of the moment of inertia of DCl being nearly twice 
that of HCl, near-resonant collisions are always possible for an HCl molecule interacting 
with a DCl molecule, whereas for DCl interacting with HCl, near-resonant collisions are 
possible for only half the energy levels. On this basis, the following values for the reso- 
nant correction were assumed: 


®HC1-DC1 " ®HC1-HC1 


(9a) 


and 


_ ®DC1-DC1 
°DC1-HC1 “ o 


(9b) 


The results of the mixture calculations are presented in table DC and are shown 
graphically in figure 2. The method referred to as ’’inelastic” uses the equation of 
Monchick, Pereira, and Mason (ref. 5). The Hirschfelder-Eucken calculations, which 
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Figure 2. - Deviation of thermal conductivities of equimolar hydrogen 
chloride - deuterium chloride mixtures from linear average. 


involve no inelastic effects, are also given in order to show the magnitude of the first- 
order correction terms. In addition, both these calculations were carried out with the 
assumption that there was no resonant exchange during unlike collisions 
^^HCl-DCl “ ®DC1-HC1 “ wovild be the case if the rotational exchange process 

actually required an exact resonance. The experimentally determined conductivities are 
also listed in table IX. The experimental points plotted in figure 2 have an apparent 
smoothness of about five parts in ten thousand. 

It is quite clear that the conductivities obtained when no vinlike interaction is as- 
sumed are much too large; however, both methods give very satisfactory agreement 
with experiment when unlike interaction is taken into accomit. In fact, the conductivities 
obtained from the inelastic calculations agree with experiment to within one part per 
thousand. Monchick, Pereira, and Mason (ref. 5) concluded that, for most purposes, 
the Hirschfelder-Eucken mixture formula is satisfactory for nonpolar gases or a 
mixture of a nonpolar and a polar gas provided that inelastic effects are included in cal- 
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dilations for the pure components. It appears that their conclusion also holds for a mix- 

3 

ture of polar gases, at least for an equimolar HCl-DCl mixture. 

The final system to be discussed is that of ammonia and its isotopically substituted 
counterparts. The thermal conductivity of equimolar NHg-NDg was measured at 300° K 
(see table in). The mixture conductivity was suprisingly low; indeed, the conductivities 
of both pure NHg and pure NDg are greater than that of the mixture. These gases ex- 
change rapidly in a glass vessel even at room temperature (ref. 43), which results in a 
four -component system: NHg, NDg, NHgD, and NHDg. Part of the discrepancy in the 
mixture conductivity can be rationalized on the basis of the lower heat capacities of 
NHgD and NHDg. 

Because of the complexity of the NHq-ND, mixture, it was decided to test the exact 

0 0^5 

resonance requirement using the equimolar NHg-N Hg mixture instead. For this mix- 
ture, the composition will remain constant despite any exchange reactions that might 
take place. From table HI, it is apparent that the mixture conductivity lies midway be- 
tween that of the pure components. In fact, the weighted average of the conductivities 
of the two components agrees with the mixture conductivity to less than one part per 
thousand. In the light of the HCl-DCl results, it seems quite clear that despite the lack 
of exact rotational energy level differences (the moments of inertia differ by about 0. 3 
percent) rotational exchange is occurring. Otherwise, the conductivity of the mixture 
would be considerably larger (about 4 percent) than the average of the pure components. 


CONCLUDING REMARKS 

The experimental results are consistent with the notion that a resonant, or near- 
resonant exchange of rotational quanta is important in deter ming the thermal conduc- 
tivity of polar gases. This is substantiated by measurements on the pure compounds and 
also on HCl-DCl and N^^Hg-N^^Hg mixtures. 


^It should be mentioned that the HCl-DCl system is actually a four-component system 

35 37 

due to the existence of two chlorine isotopes. Cl and Cl . However, the moments of 

inertia of HCl^^ and HCl^"^, as well as DCl^^ and DCl^"^, differ by only 2 to 3 parts per 

thousand. Since rotational exchange appears to be possible between HCl-DCl, there is 

35 37 

no doubt that it would occur between HCl -HCl , etc. 


13 


Results on NHg-NDg and also HgO-DgO mixtures (ref. 2) are complicated since ex- 
change reactions occur (ref. 4), and these are in fact multicomponent systems. None- 
theless, it seems likely that these data can also be explained in terms of considerable 
near-resonant exchange between unlike molecules and somewhat larger collision cross 
sections for unsymmetric species such as NH 2 D, NHDg, and HDO, 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 11, 1965. 
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TABLE I. - EXPERIMENTAL ELECTROMOTIVE FORCE 


MEASUREMENTS RELATIVE TO ARGON 


Gas 


Temperature, 

°K 



299.8 

333.0 

375. 5 

423.3 

472.3 



Electromotive force, mV 


Helium 

52. 253 

49. 342 

46. 068 

42. 822 

40. 364 

Nitrogen 

17. 700 

16. 671 

15. 552 

14. 359 

13. 583 

Methane 

28. 009 

27. 657 

27. 392 

27. 003 

26. 692 

Methane -d^ 

28. 878 

28. 941 

28.911 

28. 550 

28. 154 

Methane mixture 

28.483 

28. 353 

28.215 

27. 840 

27. 491 

( 0 . 52 CH 4 - 0 . 48 CD^) 








Temperature, 




300.0 

329. 1 

374 . 6 

424.2 

474. 5 



Electromotive force, mV 


Helium 

55.796 

53. 174 

46. 331 

43. 097 

40. 313 

Nitrogen 

18. 872 

17. 968 

15.611 

14. 509 

13. 523 

Ammonia 

18. 437 

19. 548 

19. 607 

20. 584 

21. 173 

Ammonia-dg 
Ammonia mixture 

18. 286 
17. 857 

19.755 

20. 119 

21. 276 

21.938 



19.921 


(0. 5 NHo-0. 5 ND J 
Ammonia-N 

17. 398 




Ammonia mixture 

17. 963 

20. 295 




(0. 52 NH 3 - 0 . 48 N^^Hg) 








Temperature, 




300. 1 

328. 5 

374.8 

423. 1 

471.4 


Electromotive force, mV 

Helium 

52. 375 

49. 615 

46. 040 

43. 076 

40. 504 

Nitrogen 

17.736 

16. 765 

15. 509 

14.476 

13. 563 

Hydrogen chloride 

-12. 647 

- 11.010 

- 8 . 574 

- 6 . 304 

-4. 394 

Deuterium chloride 

-12.968 

- 11 . 162 

- 8 . 564 

- 6 . 066 

-3.963 

Hydrogen chloride mixture 

- 12. 473 

-10. 749 

- 8 . 339 

- 6 . 008 

-4. 066 

(0. 5 HCl-0. 5 DCl) 








TABLE n. - THERMAL CONDUCTIVITIES OF HELIUM 
AND ARGON ASSUMED FOR CALIBRATION PURPOSES 


Temperature, 

°K 

Thermal conductivity, /ical/(cm)(sec)(^) 

Helium 

Argon 

299.8 

370.7 

42. 39 

300.0 

370.9 

42. 42 

300. 1 

371.0 

42.43 

328. 5 

394.5 

45. 50 

329. 1 

395.0 

45. 55 

333.0 

398.2 

45. 99 

374.6 

432.0 

50. 50 

374.8 

432.2 

50. 53 

375. 5 

432.8 

50. 60 

423. 1 

471.4 

55. 52 

423.3 

471.6 

55. 54 

424.2 

472.3 

55. 62 

471.4 

508.9 

60. 15 

472.3 

509. 5 

60. 23 

474. 5 

511.2 

60. 44 



TABLE m, 


EXPERIMENTAL THERMAL CONDUCTIVIES 


Gas 


Temperature 




299.8 

333.0 

375. 5 

423.3 

472.3 


Thermal conductivity, jLtcal/(cm)(sec)(^K) 

Methane 

80. 70 

91.21 

106. 55 

125. 17 

144. 48 

Methane -d^ 

83. 03 

95. 57 

113. 51 

134. 86 

156. 47 

Methane mixture 

81.96 

93. 52 

110.21 

130. 23 

150. 80 

(0. 52 CH4-0. 48 CD4) 
Nitrogen 

60. 56q 

65. 59 

72. 10 

78.87 

85.65 



Temperature 




300.0 

329. 1 

374.6 

424.2 

474. 5 


Thermal conductivity, /ical/(cm)(sec)(^K) 

Ammonia 

59.96 

67. 51 

80.64 

96. 13 

112. 58 

Ammonia-dg 
Ammonia mixture 

59.76 
59. 19 

67. 86 

81.91 

98. 54 

116. 20 



93. 93 


(0. 5 NH«-0. 5 NDo) 
Ammonia-N 

58. 60 




Ammonia mixture 

59. 33 

95. 16 




(0. 52 NHg-0. 48 N^®Hg) 






Nitrogen 

60. 55^ 

64.98 

71. 89 

79. 12 

85. 83 



Temperature 




300. 1 

328. 5 

374.8 

423. 1 

471.4 


Thermal conductivity, jLtcal/(cm)(sec)(^K) 

Hydrogen chloride 

34.95 

38. 03 

43. 39 

49. 17 

54.90 

Deuterium chloride 

34. 80 

37.95 

43. 40 

49. 38 

55. 37 

Hydrogen chloride mixture 

35. 04 

38. 18 

43. 56 

49. 34 

55. 26 

(0. 5 HCl-0. 5 DCl) 
Nitrogen 

60. 6O2 

64. 90 

71.92 

78.91 

85. 35 



TABLE IV. - MOLECULAR PROPERTIES OF TEST GASES 


Gas 

Moments of inertia, 
g-cm^xlO^^ 


Ia 

% = ic 

Methane 

^5. 3326 

= ^A 

Methane -d^ 

‘’lO. 7073 


Ammonia 

4140 

2. 8087 

1 5 

Ammonia- N 

‘^4.4140 

2. 8159 

Ammonia-dq 

‘^8. 7960 

5. 4129 

o 

Hydrogen chloride 


®2. 6431 
^5. 1404 

Deuterium chloride 



Molecular 

weight 


16, 0425 
20. 0678 
17. 0320 
18. 0330 
20. 0500 
^36.4693 
^37. 4763 


Dipole moment, 
18 

esu-cmxio 


0 

0 

*^1.477 
®1.477„ 
^1. 509 
^1.081 
^1. 085 


^Ref. 26. 

*’Ref, 27. 

‘^Calculated using computer program described in ref. 44. 


‘^Ref. 45. 

^Interpolated from data for NHg and NDg on basis of the reduced 
mass. 

W. 29. gg 

^Weighted average of the molecular properties of HCl and 


*^Ref. 46. 

3 5 

^Weighted average of the molecular properties of DCl and 
DCl^"^. 



TABLE V. - DATA ANALYSIS FOR PURE METHANE AND METHANE-d^ 





Temperature, 





299.8 

333.0 

375.5 

423.3 

472.3 

Viscosity of methane, /xP 

112.5 

123.3 

136.0 

150.0 

163.6 

Ci„t/R 

Methane 


1.80134 

1.96780 

2.22166 

2. 54276 

2.89179 

Methane -d^ 

1 

2.39526 

2.70998 

3. 12661 

3. 59446 

4.06090 

pDjj/r,(6/5(A*» 
Collision number, Z . 

1.3190 

1.3187 

1.3184 

1.3194 

1.3207 

Methane 


4.0 

3.6 

3.8 

3.6 

3.0 

Methane -d^ 


3.4 

2.8 

2.7 

2.4 

2.0 

^rot, (CD 4 ) 


0.86 

0.79 

0.71 

0.67 

0.66 

^rot, (CH 4 ) 
^rot, (CD 4 ) 


0.65 

0.65 

0.65 

0.65 

0.65 

^rot, (CH 4 ) 

cal, eq. (6) 


TABLE VI. - DATA ANALYSIS FOR PURE HYDROGEN CHLORIDE 
AND DEUTERIUM CHLORIDE 





Temperature, 





300. 1 

328. 5 

374.8 

423. 1 

471.4 

Viscosity of hydrogen chloride, p.P 

146.4 

160.7 

183.8 

207.3 

229.0 

Cint/« 







Hydrogen chloride 

1 . 00430 

1.00501 

1 . 00704 

1.01116 

1.01834 

Deuterium chloride 

1.00859 

1.01346 

1.02658 

1. 04799 

1. 07709 

PDjj/t7(=6/5(A*)) 

1. 3254 

1. 3243 

1. 3232 

1.3214 

1.3190 

Collision number, 






Hydrogen chloride 

6.2 

4.6 

3.6 

3.2 

3.0 

Deuterium chloride 

2.6 

2.3 

2.0 

1.9 

1.9 

^rot, (DO) 


0.43 

0. 50 

0. 56 

0.60 

0. 62 

^rot, (Ha) 







^rot, (DCI) 


0. 54 

0. 54 

0. 54 

0. 54 

0. 54 

5rot, (HC1)_ 

cal, eq. (6) 






Resonant correction, 6 






Hydrogen chloride 

0.48 

0, 42 

0.34 

0.28 

0, 24 

Deuterium chloride 

0.18 

0.16 

0.13 



0.10 

0.09 



TABLE Vn. - DATA ANALYSIS FOR PURE AMMONIA 


AMMONU-N^^, AND AMMONIA -dg 






Temperature, 





300.0 

329.1 

374.6 

424.2 

474. 5 

Viscosity of ammonia, iiV 

101.53 

112.2 

129.0 

147.2 

165.8 


Ammonia 


1.78007 

1.87133 

2.02926 

2.21460 

2.40886 


Ammonia - 
Ammonia- 

n15 

2. 10154 
1.78205 

2.25252 

2.49342 

2.75482 

2.21792 

3.01423 






pD /t7(=6/5(A*)) 

1.3202 

1.3218 

1.3241 

1.3254 

1.3243 

Collision number, Z . 







Ammonia 


2.3 

2.1 

1.9 

1.7 

1.6 


Ammonia- 

^3 

1.6 

1.5 

1.3 

1.3 

1.2 


15 

Ammonia -N 

2. 5 



1.8 







r 

^rot, (NDg) 


0.69 

0.70 

0.72 

0.74 

0.77 


^rot, (NHg) 



^rot, (NDg) 


0.64 

0.64 

0.64 

0.64 

0.64 


^rot, (NHg) 

cal, eq. (6) 

Resonant correction, 6 







Ammonia 


0.15 

0.12 

0.09 

0.07 

0.06 


Ammonia -dg 
15 

Ammonia -N 

0. 05 
0. 15 

0.04 

0.03 

0.02 

0.07 

0.02 






TABLE Vm. - COMPARISON OF EXPERIMENTAL METHANE - METHANE-d^ 
MIXTURE THERMAL CONDUCTIVITIES WITH THOSE 
CALCULATED FROM LINEAR MIXING RULE 




Temperature, 




299.8 

333.0 

375.5 

422.3 

472.3 


Thermal conductivity, jLi cal/(cm) (sec) (°K) 

Experimental (0. 52 CH^-0. 48 CD^) 

81. 96 

93. 52 

110.21 

130.23 

150. 80 

Calculated mixing rule 

p 

00 

00 

93.32 

109. 91 

129.85 

150. 27 

(0. 52 CH^-O. 48 CD^) 

Ratio of experimental to calculated, 

1.0017 

1.0022 

1. 0027 

1.0029 

1.0035 

^^exp/^cal^ 







^he pure gas conductivities are given in table HI. 



TABLE rx. - COMPARISON OF EXPERIMENTAL AND CALCULATED 
VALUES FOR THERMAL CONDUCTiyiTY OF EQUIMOLAR 
HYDROGEN CHLORIDE - DEUTERIUM CHLORIDE MIXTURES 


Method 


Temperature, 



300. 1 

328. 5 

374.8 

423. 1 

471.4 


Thermal conductivity 

jLtcal/(cm)(sec)(^K) 

. 


Inelastic, unlike interaction 

35.08 

38.21 

43. 62 

49.49 

55. 33 

Inelastic, no unlike interaction 

36.35 

39. 59 

45.09 

50. 99 

56. 80 

Hirschfelder-Eucken, unlike 

34.98 

38. 09 

43.48 

49.36 

55. 22 

interaction 






HirschfeLder-Eucken, no unlike 

35. 58 

38. 64 

44.01 

49. 88 

55.75 

interaction 






Linear mixing rule 

34. 87 

37. 99 

43. 39 

49.28 

55. 13 

AexpHci ^ ^expj3Cl\ 

\ 2 / 






Experimental 1 

35. 04 

38. 18 

43. 56 1 

49.43 

55. 26 


NASA-Langley, 1966 E-2883 
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